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Abstract 

Arguments are presented to the effect that embedding semi-ffexibie 
(wormlike) ideal polymers into a fluctuating, flexible surface leads 
to an effective attractive orientational interaction between polymer 
segments that precipitates an orientational ordering transition of the 
polymer chains on the surface even in the case of otherwise ideal (non- 
interacting) chains. The orientational interaction is analogous to the 
(zero order) Casimir force and is due to the energy change in surface 
conformational fluctuations in the presence of embedded semi-flexible 
chains. February 1, 2008 
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1 Introduction 



Polymers or polymer-like structures confined to lye in an embedding surface 
can appear in the context of different membrane and/or interfacial phenom- 
ena. Partially polymerized membranes of unsaturated amphiphiles can some- 
times make long trains of joined monomers that behave as surface - embedded 
polymer chains M. Hydrophobic polymers in aqueous solutions can become 
embedded or multiply attached to hydrocarbon interior of membranes Jj]], 
while certain polymers can also form monomolecular films on the water-air 
interface fl2|. It has been observed that embedded polymers in the case of 
partially polymerized membranes can promote wrinkling transition in mem- 
branes || while multiply anchored polymers give rise to bulging and budding 
of membranes 0] . In the case of interface polymer films || it has been shown 
that they make segregated domains even at very low (submonolayer) concen- 
trations where one would expect the polymers to be evenly distributed over 
the surface. 



This plethora of new phenomena that appear closely connected with 
the fact that polymers are not only adsorbed but actually embedded into a 
supporting surface have recently become a focus of several theoretical studies 
U . In this contribution we try to find out what is the influence of finite tem- 
perature, which in general leads to shape fluctuations in flexible (membrane) 
surfaces or interfaces, on the properties of embedded polymer chains. In the 
framework of a semiflexible polymer chain model we show that finite tem- 
perature leads to membrane fluctuation generated effective orienting forces 
between different polymer segments. This interaction formally resembles the 
usual nematic interaction used in various models of the nematic ordering of 
long flexible chains fl2| and is thus not surprising that it leads to an orien- 
tational ordering transition. 



The interesting point is that this transition is driven by shape fluc- 
tuations of the supporting surface into which the polymer is embedded. 

The problem we address here is related to the recent investigations 
of fluctuation-induced forces between manifolds immersed in a correlated 
fluid |J. The existence of surface fluctuations promoted attraction between 
different segments of the embedded polymer chain can also be related to the 
zero - order Casimir forces (Keesom forces) between beads attached to flexible 
membranes M. Other related problems include adsorption of polymers on 
soft, fluctuating surfaces M or the collapse of polymers on fluctuating surfaces 



The outline of the paper is as follows: I shall first derive the partition 
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function of a surface embedded semiflexible polymer chain at zero temper- 
ature. Then I shall relax the zero temperature constraint and derive the 
effective partition function for a polymer, after the supporting surface de- 
grees of freedom have been integrated out. At last an approximate solution 
for this effective partition function will be given and discussed. 



2 Analysis 



2.1 Semiflexible polymer chain embedded in a surface 

I use a semiflexible (worm - like chain) model to describe the polymer chain. 
Though the theory presented below is formulated for a single long polymer 
chain, it can be trivially generalized to a many - chain system by assuming 
that the variable n has a continuous part N a , describing a single chain, and 
a discrete part describing different chains a. 

I investigate the statistical properties of a semiflexible polymer chain 
embedded in a fluctuating surface. Let us assume that the n-th bead of 
the polymer has coordinates r(n) = (x(n),y(n), z{n)) and that the quenched 
profile of the embedding surface can be written in a Monge parametrization 
z = y) = C(p), where I introduced the two dimensional vector p = (x, y). 
The embedding ansatz for the polymer can thus be written in a form 

*(n) = C(p(n)). (1) 

I also assume that the polymers can be described in terms of the semiflexible 
chain model, where the configurational energy is given as 

rN ( d?r(n)\ 2 r N 
0W„(r(n)) = §/fe jf ( dn = \(3e jf (i{n)f dn (2) 

where /3e = -p, with I being the persistence length and N is the number 
of monomers. For a differentiable curve the tangent r(n) is a unit vector, 
r(n) 2 = 1. The partition function of a semiflexible chain embedded in the 
surface (,(x,y) can thus be written in a succinct form 

5 (C (p(n))) =/.../ Vv(n) J[ 6 (z(n) - ((p(n))) ]J 6 (r(n) 2 - l) ^ H ^\ 

(3) 

The explicit incorporation of the constraint that the local tangent has to be 
a unit vector permits us to perform the functional integration in Eq|J over 
an unconstrained set of r(n). 
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By integrating out the z(n) variables and introducing an integral 
representation for the delta function I obtain the partition function in the 
form 

H (C (p(n))) = /•■•/ Vp{n)VX{n) e ~ W(PM,C(p(«))) (4) 
where the effective configurational Hamiltonian reads 

rN rN 

/9H(p(n), ((p(n))) = |/3e jf (p(n) 2 + C(pW) 2 ) dn + iJ q A(n) (p(n) 2 + C(pW) 2 ) 

(5) 

It is at this point where a first major approximation is necessary in order to 
proceed. Instead of enforcing the condition r(n) 2 = 1 locally at each point 
along the polymer chain, we presume that it is only satisfied globally [[HJ i.e. 
on the average 



;r(n) 2 ) = (p(n) 2 + C(p(n)) 2 ) = l, (6) 

where the average is over all polymer configurations. Thus the magnitude of 
the tangent can fluctuate locally but on the average, over the entire polymer 
chain, it is fixed. Replacement of the local continuity constraint with a global 
one Eq. |6| can also expressed as A(n) — ► A = const.. Thus one remains (after 
making at the same time also the transformation zA — > A) with the following 
form of the partition function 

5 (C (p(n))) = J...J Vp{n) exp (-/3W(p(n))) , (7) 

where some irrelevant constant terms have been omitted. Also 

rN 



ply 

[3H{p{n)) = \(3e J (p(n) 2 + C(n) 2 ) dn + 

rN 

+ \j o {p(nf + C(nf) dn 



(8) 



where we have limited ourselves to the terms of the second order in the 
derivatives of p(n). 



2.2 Conformational fluctuations of the embedding sur- 
face 



Going now to a non-zero temperature one has to take into account that the 
embedding surface is allowed to fluctuate, taking the polymer with it on its 
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conformational wanderings. The unrestrained partition function in this case 
is 

5 (AO = J...J Vp(n)V( (p) exp (-(3Ht (C (p) , p(n))) (9) 
where after introducing the Fourier transform £(p) = J2q C(Q) e iC *P 

I3H T (C (p) , p{n)) = \pe p{nfdn + X p{nfdn + 

Jo Jo 

IEC(QMQ,Q')C(-QO + ^EnC(Q)) 

Q,Q' Q 

(10) 

where 

r N 

P(Q, Q'; p(n)) = (3e dn (Qp{n)) (Q'p(n)) e^-WW + 

JO 

/-iV 

+ 2A / dn (Qp(n)) (Q'p(ra)) e <Q-Q')PM 
Jo 

V(£(Q)) is the intrinsic deformation energy of the surface describing the 
energetics of the surface in the absence of embedded polymers, I assume it 
to be of the form 

V(C(Q)) = V(Q) |C(Q)| 2 - (12) 
With this ansatz one can integrate out the ((p) degrees of freedom obtaining 

S(p(n)) = /. . ./l>C (P) S(C (p),p(n)) = 

r f rN rN N 

.y Pp(n) x exp(-i/3e jf p(n) 2 dn - XJ p{nfdn - \Tr In £(Q, Q'; pin)) 



(13) 



with 

0(Q, Q'; P(n)) = (pV(Q)) 5(Q - Q') + £>(Q, Q'; p(n)). (14) 

At this point I use the standard identity relating Tr In of an operator with 
a coupling integral of its resolvent 

Tr ln£(Q,Q';p(n))=Tr [' dfiK^Q, Q'; p(n)), (15) 

jo 

where the resolvent 7^(Q, Q'; p(n)) is defined as 



^(Q,Q';p(n))=P(Q,Q , ;p(n)) x (j + /iX>(Q, Q'; p(n)) (^(g))" 1 

(16) 

with X being the identity operator in Q space. I have omitted a term that does 
not depend on the coordinates p(n) from Eq.|15]. Evaluating the resolvent to 
the first order in [i I obtain 

Tr ln£(Q,Q^p(n)) = TrX?(Q,Q';p(n)) {(5V{Q)y x -\ Tr (Z>(Q, Q'; p{n)) ((3V (Q)y 1 ) 2 + . . . . 

(17) 
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What one remains with at the end is an effective interaction between polymer 
segments as if they would be confined to a planar surface. It has the form 



\Tr lnC?(Q,Q';p(n)) = ^ [ dnp t (n)p k (n) £ + 

f N Q Q 



Q 



( 2A ) 2 f N f N dndn'b(n)b,(n)b,(n')b (ri) V e ^P( n )-P^')) y -»q'(p(n)-p(n')) , 

4 io Jo P,WP*WPH" jp m (n j ^ e ^ (/3V(Q')) 

(18) 

I have limited myself to second order terms in the derivatives of p(n) at any 
n (self - energy term) and to fourth order terms for two different values of n 
(interaction term). 

Introducing now the orientational correlation function between the 
directions of the surface normals n(p) in the absence of polymer chains and 
for small surface deformations as 

<n,(pK(p')> = ( d -^ d 4f-) = ~ E Q l Q' fe (C(Q)C(Q0) e^'^ 

\ U Pi U Pk I Q,Q' 

(19) 

where due to the harmonic ansatz for the configurational energy of the bare 
surface I have 

kT 

viQ) 

thus 



(C(Q)C(Q')) = ^T777^ 5 (Q + Q'), (20) 



(3 (n,(p)n fe (p')) = \m*{p ~ P') = \ E e^P P'\ (21) 

The final expression for the partition function of a surface - embedded poly- 
mer chain at finite temperature is thus 

E (iV) =/•••/ Vp{n) e -Pn eff (p(n),p(n)) (22) 

with 

rN rN 



rl\ rl\ 

/3H e ff(p(n),p(n)) = \j3e p(n) 2 dn + X p(n) 2 dn + 

rN 

+ f3e / dnp^nYp^n) (n,(0)n fc (0)) + 
jo 

+ 2X f dn Pl (n)p k (n) (n,(0)n fc (0)> - 
Jo 

rN rN 

-AX 2 / dndn 1 ' p^n) p k {n) p^n') p m (n') (n i (n)n m (n')) (n fc (n')n,(w)> + . . . 
Jo Jo 



(23) 
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Obviously the integration of the surface degrees of freedom leads to an effec- 
tive attraction between the polymer segments, which is orientation depen- 
dent. This attraction has the same origin as the usual (zero order) Casimir 
force [n and we could simply dubb it the nematic Casimir force, or maybe 
more appropriately the nematic Keesom force. The fluctuations of the em- 
bedding surface play in this case the role of vacuum fluctuations of the elec- 
tromagnetic field 



2.3 Effective embedded - polymer partition function 
and orientational ordering 



Introducing now the orientational tensor a ik (p) as J 1 

rN 



<T ik {p)= dn p^njp^n) 5{p - p{n)) (24) 
Jo 



I can write the partition function alternatively as 

rN 
10 



E (iV) = J...J Vp(n)Va lk (p) <*( jf dn Pi (n)p k (n) 5(p - p{n)) - a tk (p)j e HW.,,(P(»),P(nW(P)) , 

(25) 



where 



/3H e ff(p(n),p(n),a ik (p)) = \(3e p(nfdn + Ay p{nfdn + 

rN 

+ \(3e J dn ^(n)p fe (n) ^(0) + 
+ \Jd 2 p a ik (p) ^(0) - 

- A 2 J J d 2 pd 2 p' a ik (p)a lm (p') T im {p - p')T kl {p - pi) + . . . 

(26) 

In order to proceed from here one has to specify the form of V(Q), valid 
for the bare surface. To get some general trends and estimates I shall limit 
myself to the following simple form of !Fi k (p) 

T ik {p) = 5 ik F{\p\). (27) 

Also writing the delta function in an integral representation I obtain the 
following form of the effective Hamiltonian 

rN 

/3n(p(n),p(n),a ik (p),tJ; ik (p)) = (3n e ff(p(n),p(n),a ik (p)) + i\ dn p^njp^n) ip ik (p(n)) 
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- % J d 2 p ipik(p)<?ik(p) = 

rN rN 

= ±/3e(l + F(0)) / p{n) 2 dn + X p{n) 2 dn + 
Jo Jo 

+ (AF(O)) J d 2 p Tr{a lk (p)) - 

- X 2 J J d 2 pd 2 p' Tr(a lk (p)a kl (p'))F 2 (\p-p'\) + 
dn p l (n)p fe (n) ip ik (p(n)) - i J d 2 p i) ik {p)(r ik (p) 



+ i 



and 

3 (N) =[■ ■ ■[ T>p(n)Va ik (p)Vip ik (p) exp (-/3H(p(n), p(n),a ik (p), ip ik (p) 

(29) 

The above partition function can not be evaluated explicitely because i/ji k (p(n)) 
depends on the polymer coordinates and the functional integrals with respect 
to p{n) and ipi k (p(n)) can not be separated. However, as recognized by Gupta 
and Edwards [I2[ another approximation suggests itself. Replacing the fields 
Cik(p) and 4>i k (p) by their average values a ik and ipi k , independent of p ( and 
at the same time making the transformation %^i k — > ipi k ), one can first of 
all explicitely evaluate the p(n) part of the functional integral 

j T>p(n) e"^ £{1+F(0)) J N P(n) 2 dn - \^p(nfdn - ^ dn p t (n)p k (n) fa = 



(30) 

where 

M lk (j)= (|/5e(l + F(0))(^) 4 + A(^) 2 )^ fc + (^) 2 ^ fc . (31) 
Above I introduced the Rouse modes of the polymer chain as 

j=+oo 

p(n)= £ p(j)e 2 ^. (32) 

j=-oo 

The sum of In DetAii k (j) can be evaluated explicitely and in the limit of long 
chains where the sum over j can be replaced with an appropriate integration 



one obtains 12 



J= + OD 

\ E lnDetM ik (j) = f E 

j=-oo 



+ x 



\ f (i + JW 



(33) 



(28 
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where a is the index of the eigenvalue of ipik- I have omitted a divergent 
contribution to the above expression since it is irrelevant for subsequent de- 
velopements. In the limit of 'mean fields' one can now write the free energy 
corresponding to the partition function Eq.p9] as 



MF 



f £ 



A 



*\ f (i + no)) 



-NX-SJ2 *«$a + AF(0)5 °a - A 2 5a 2 £ a 2 a +. 

a a a 

(34) 



where S = J d 2 p is the projected area of the surface and I introduced the 
second virial coefficient a 2 as 



a 2 = Jd 2 pF 2 (\p\). 



(35) 



Since the mean - field free energy is a function of (5Tmf = 0J~mf(X, &a, ^«) 
I get the equilibrium solution by minimizing it with respect to all three 
variables: 

1 



MF) 



8(J3Fmf) 
da a 

diP^MF) 
dX 










cr, 



N 

s 



4^f(l + F(0)) v^TA 
4 = AF(0) - 2A 2 a 2 a a + 0(a 2 a ) 

« 1 AT 



4Jf(l + F(0)) vC+A 5 



F(0) ]>> a + 0(a 2 



(36) 



In order to be consistent the above equations should be linear in the field a a . 



These equations look very much like the analogous equations ob- 



tained by Gupta and Edwards |L2[ for athermal flexible polymer chains with 
orientation dependent attraction. Solving these equations explicitely for the 
three fields I obtain the following results. First of all a a is a solution of 



[a x — a,. 



2a 2 X 2 - 



(« 



v- 



(a x + a v 



\ 



8((3e) (1 + F(0)) ola* 



From the last equation I obtain 

? CTq = (1 + F(0)) 
which is simply a consequence of the continuity of the chain since 



N 

s 



(37) 



(38) 



Tr a a 



N 

~S 

N 



P [Ji, 



N 
~S 



(i - (i 2 H)) 



'dC(p) d((p). 



dpi dp h 



Pi(n)pk(n) 



N 
~S 



1 - F(0) (p 2 {n) 



(39) 
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These are the two equations that have to be solved consistently in order to 
get the dependence of the eigenvalues of the orientation tensor on the surface 
monomer density, <j. There is also a boundary condition that we have to 
take into account, viz. that in the case of a stiff surface F — > the polymer 
statistics should reduce to the case of a semiflexible chain embedded in a 2-D 
surface. 



The appropriate solution for A can now be derived from the first two 



equations of Eqs.BBI in the form 
(1 + F(0)) 2 



A 



N 
S 



V 



(gg) (I 
\r 4 e (i+F(o)r 



E- 



(40) 



The nonisotropic (i.e. a x 7^ a y ) solution for the orientational tensor can be 
obtained from 



S 



(a x + a y ) 16(/3e)(l + F(0))a 2 A 2 
Defining the orientational order parameter 5 as 



a 



[1±S) 



(41) 



(42) 



;i + f(o)) 

so that Tr a a satisfies the constraint Eq. |38|, one can derive an approximate 
equation for 5 that has two solutions. One corresponding to the isotropic 
state (5 = 0) and the other one corresponding to the orientationally ordered 
phase of polymers (S 7^ 0), 







S 



\ 



. (i+^(Q)) 2 

' (/3e)A 2 02 f 

(i+^(o)) 2 . (i+p(Q)) 2 

(/3e)A 2 a 2 f 



(/3e)A 2 a 2 



.f 



> 1 



< 1 



(43) 



The above approximate solution is valid only in the vicinity of the critical 
point 5 = 0, where the general solution for A in terms of 5 can be approx- 
imated with its lowest order term. Away from the critical point additional 
terms in the expansion of A = A(5) would have to be taken into account. 
The discussion in this contribution will not be extended to this regime of 
orientational order parameter. 



3 Discussion 



As already stated the existence of effective attraction between segments of a 
polymer chain embedded in a flexible surface at non-zero temperature is not 
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completely unexpected. D'Hoker et al. || have found an effective (Casimir) 
attraction between beads embedded in strings and membranes. Though their 
model applies to a different situation than our own calculation, its conse- 
quences stemm from the same basic physical mechanisms. 

The effective nematic interaction of the form Jq t J^ dndn' ' p^ri) p k (n) 
Pi(n')p m (n') (n i (n)n m (n / )) (n fc (n')n^(n)) can in fact be viewed as a special 
type of Casimir force between different segments of a polymer chain. It is due 
to the change in the energy of conformational fluctuations of the supporting 
surface ( 'zero point energy' in the case of electromagnetic Casimir forces) 



in the presence of embedded chains. It is analogous to the zero order van der 
Waals interaction energy that is due purely to thermodynamic fluctuations 

0- 

The dependence of this Casimir force on the orientation of the poly- 
mer segments comes essentially from the geometric constraint for the polymer 
chain, i.e. r 2 (n) = 1, since some of the 3-D configurations are excluded if the 
chain is confined to lye in the embedding 2-D surface. Its spatial dependence 
on the other hand, thus the form of F(\p\), stemms purely from the nature 
of the bare surface fluctuations. 



Once the presence of a nematic Casimir interaction is derived it nat- 
urally leads to an orientational ordering transition. Also the order of the 
transition clearly follows from the fact that the effective polymer partition 
function corresponds to a polymer confined to a 2-D surface. Nematic tran- 
sition on a 2-D surface is of a second order [HJ. Whether the order of the 
transition is preserved for non-ideal polymers and/or higher orders in the 
perturbation expansion remains to be assesed. 



The present calculation is of importance for phenomena involving in- 
terface embedded polymers. Monomolecular films of poly(dimethyl siloxane) 
H show domains of differing surface density even below the submonolayer 
concentrations. Orientational ordering described in this contribution could 
possibly lead to domains of different surface density of polymers at quite 
small average surface densitites. One way of testing this mechanism is to 
quench the free surface fluctuations and thus attenuate the nematic Casimir 
attraction which should prevent the ordering of the embedded polymers. 

Some of the open questions that we leave for further discussion is 
whether the nature of the transition is conserved also in higher orders of the 
perturbation expansion in terms of F and what is the interplay between the 
repulsive (steric) interactions and the nematic Casimir attraction. 
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